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Abstract 

A simple and innocent modification of Poisson's equation leads to a 
modified Newtonnian theory of gravitation where a localized and positive 
energy density of the gravitational field contributes to its own source. The 
result is that the total active gravitational mass of a compact object is the 
sum of its proper mass and an evanescent gravitational mass which is a 
mass equivalent to the gravitational energy. 

Conventional wisdom subordinates the concept of the energy of the Newton- 
nian gravitational field to the potential energy of the source. As if the unique 
interest of the first was to provide us with an alternate method to calculate the 
second. We are thus led to accept that the energy of the gravitational field is 
negative. 

Since the potential function V is defined only up to an additive constant, 
conventional wisdom has also made us comfortable with choosing this constant 
such that Vqo is zero, while identifying it with the square of a velocity opens 
new opportunities of interpretation. 

The first section below introduces the notations, develops somewhat the 
preceding remarks and serves as a transition to the second section. There we 
define the energy density ag of the gravitational field as a localized positive 
quantity, unrelated a priori to the potential energy density, and modify the 
Poisson equation The new equation is not linear and is formally equivalent to 
the Poisson equation with a source term that is the sum of the proper mass 
density pp and the evanescent gravitational mass one pg = ag/c^, where Cg, is a 
new parameter with the physical dimensions of a velocity, that may be related 
to Voo or not. 

An obvious consequence of this modified theory is that the active gravita- 
tional mass is always greater than its proper mass. The gain factor might have 
any value greater than 1, depending on the object under consideration. 

The last section contents a few elements of comparison between the plain 
Newtonnian theory and the modified one, as well as some comments about 
possible choices of the new parameter Cg. Its value should be derived from 
an appropriate confrontation of the new theory with observational data, an 
undertaking that is well beyond the scope of this paper and the competence of 
its author. 

Newtonnian gravity.- Newtonnian gravity is based on Poisson's equation: 
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Ay = 47rGpp (1) 

This equation being Hnear the gravitational potential V is defined only up to 
an additive constant that it is usually chosen so that V tends to zero at infinity. 
We shall see in a moment that making a different choice opens interesting new 
perspectives. 

We shall restrict the scope of this paper to the case where the source of 
the field is a spherical body with finite radius R and proper mass density pp 
depending only on the distance r at the center. Then we have: 

d^V 2dV ^ ^ 

Setting: 

W^r^^ (3) 
dr 

Equ. ^ is equivalent to the couple of equations: 

— ^^nGppr\ (a) ^ = ^ (b) (4) 

whose solution is, requiring the value of W at the center to be zero: 

f f°° W 

W = 4ttG ppT^ dr, (a) V = - —dr + Voo (b) (5) 

Jo Jr ^ 

where Voo is an arbitrary constant. Let us call proper mass of the source the 
quantity Mp-. 

/>oo 

Mp = 4tt ppr"^ dr, (6) 
Jo 

and active gravitational mass the quantity Ma in the formula: 

GMa = hm W (7) 

r— >oc 

Mp and Ma are different physical quantities, but it follows from ([S^) that nu- 
merically we have Ma = Mp. 

We define the energy density of the gravitational field as the positive quan- 
tity: 

and using ([T]) and integrating by parts the integral giving the total energy we 
obtain: 

poo -t 

Eg^Air ay dr = —iVW)oo - 2^ / ppVr^ dr (9) 
Jo Jo 

If we choose to set Voo = then we have Eg — ~Up/2, where Up is the gravita- 
tional potential energy0: 



^See for instance T III. 2 



2 



Up^Air Vppr^ dr (10) 
Jo 

But otherwise we can write the preceding formula ^ as: 

Ma = -^{2E, + Up) (11) 
A second natural choice would be to set: 

V^= I —dr (12) 







in which case one has Vq = at the center and V positive all the way for r > 
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We shall assume that Voo is positive and write Voo = c^, so that the quantity: 



P.-'-i (13) 

^9 



will be a mass density equivalent to the gravitational energy density. 

One of the potential benefits from choosing Voo satisfying (fT^ is that the 
gravitational active mass-function: 



Ma{r) = 4^ / {2pg + {V/cl)pp)r' dr (14) 

JQ 

is then a positive monotonously increasing function of r tending to Ma as r goes 
to oo. 

The formula (jlip would be a potentially interesting one if we could relate 
Up to some luminosity mass. But since this is doubtful we proceed to the next 
section. 



Ponderous gravitational energy - In this section we consider the following 
modification of Eq. ([T]) : 

Ay = AttG{pp + pg) (15) 

where we maintain the definitions of Ug and pg in ([S]) and (|13p . We shall maintain 
also the definitions of Mp in ^ and Ma in ^ , while the discussion about the 
plausible values of the parameter Cg will be made latter. 
Eq. (fT5|) is equivalent to the couple of equations: 

— ^4nG{pp + Pgy, — ^- (16) 

The first one is a Ricatti equation that we integrate assuming that pp is constant 
and requiring {Wi/r'^){0) = 0. The interior solution Wi, i.e. for r < i?, is then: 



ar 



where: 



W, = 2rc'g [I + —Fir) ) , (17) 



Fir) ^ + ■^i^ + ^^^^r) ^ ^^g^ 

1 — sin ar cos Q;r(cos ar — 1 — sin ar) 
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s 



or 



and: 



F(r) = -cot(|r) 



While the exterior sohition We, i.e. for r > R, is: 

We = 



1 + 2Cer 



(19) 



(20) 



(21) 



where Ce is a constant of integration that has to be fixed requiring the continuity 
of W across the boundary of the source: 



W^{R)^We{R) 

On the other hand, using the definition ([7]) of Ma, we have: 
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so that setting: 



GMa = 7f 



GM„X 



(22) 



(23) 



(24) 



where X = Ma/Mp is the gain factor, we can calculate this quantity using (|17p . 
(PT|) and (im, obtaining after some elementary algebra: 



X 



8^/3 + 12VAF(A) 



3AVAF(A) 



where A is the compactness parameter. 



and: 



F(A) = 



_ 2GMp 



1 + cos VSX — sin V3A 



1 — cos V3A — sin VSX 
X is a monotonous increasing function of A becoming infinite for: 

These are the values of X for A = 1, 2, 3: 



(25) 



(26) 



(27) 



(28) 



A 


X 


1 


1.43 


2 


2.53 


3 


11.20 



^Courtesy of J. M. Aguirrcgabiria and J. Martin 
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For small values of A, X can be approximated to: 



Final comments - To compare the two theories, the Newtonnian and the 
modified one, we consider a test particle moving freely in the exterior gravita- 
tional field, r > i?, of the spherical body. Let oq be the acceleration according 
to Newtonnian theory and let oi be that one according to the modified theory: 

Defining the anomaly e as follows: 

ao 

we have, using the available definitions: 



(31) 



4r / 4r \ ^ 
e = — 1 + — — - 1 (32) 
XR \ XXR J ^ ' 

The larger value that this anomaly e can be at a distance r is: 

4 r 



^max 



A R ^ ^^^^ 

where Xmax is the value given in 

For small values of A, using the approximate expression ([2^ in ([5^ and 
approximating it again we get: 



_3_ IR 

10 ^ 47 



A (34) 



At this point, to push further the comparison between the two theories 
and eventually when trying to falsify the modified one will requires to face the 
problem of choosing Cg. This is something that we shall not attempt to do it 
here. But we comment below several open possibilities: 

(i) Choosing Cg = c is a natural choice nowadays. It would not have been 
in the 19th century. With this choice two bodies with the same compactness 
parameter A but quite diff'erent visible parameters Mp and R will have the same 
gain factor X and give raise to the same qualitative gravitational field. This 
is the case with th Sun, A = 4.3 10"^, and the huge galaxy NGC 4472, A = 
3.8 10~^, which has a proper mass of the order of 10^^ solar masses and a 
diameter of about 50 kpc. It looks surprising that two objects so different from 
the morphological point of view could share some essential features of their 
gravitational fields. 

R being the radius of the sun let us assume that r « 215 x R which is 
approximately 1 AU. We have then, calculating the anomaly defined in (|32p , e = 
1.2710^^. For much larger values of r the anomaly will remain approximately 
constant to e = 1.2810"^. 
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(ii) Choosing Cg to be a universal constant different from c would be a major 
cliange in our understanding of gravity, but we shall have to wait until observa- 
tional bounds on the value of Cg become available. This possibility does not cure 
the embarrassing precedent remark coming from the similitude of the values of 
A for the Sun and NGC 4472 

(iii) considering Cg as an invisible phenomenological parameter depending 
on the morphology of the object under consideration severely restricts the pre- 
dictability of the modified theory, but it may be useful to to take into account 
the morphology of the source of the field and fit scandalous anomalous data. 

(iv) Maintaining the definition (fT2|) of Voo in the modified theory and setting 
Cg = Voo leads to the equation: 

/•-R /'oo 

= / w^dr+ We dr (35) 

Jo JR 

Using (I17p and ([?!]) drops out and calculating the integrals we get: 
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^ <^l + tanW)A ^ /l + 2Ci^. ^3^^ 



-2i?2 I _^ tan^ 

1 tan2(ai?/2) ) ' V 2Cei? 
that, using the available definitions of a. A, and X, gives: 



1 >^^y 3A 

'^-) 4sin-(V3A/2) ^^ ^''^ 

which is an equation for A. Its unique acceptable solution is A = 1.126342730. 
With this choice Cg remains arbitrary and A takes a fixed value independent 
of the visible parameters. The corresponding common gain factor being: X = 
1.513581269 

(v) Choosing Cg{t) a function of time would allow to discuss possible evolu- 
tionary effects. 

(vi) Maybe the best choice: Cg = oo 
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